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Abstract 

The Lewenstein-Sanpera decomposition for a genèric two-qubit density matrix is obtained 
by using Wootters's basis. It is shown that the average concurrence of the decomposition is 
equal to the concurrence of the state. It is also shown that all the entanglement content of the 
state is concentrated in the Wootters's state associated with the largest eigenvalue Ai of the 
Hermitian matrix J~^fpp\fp ■ It is shown that a given density matrix p with corresponding set 
of positivc numbcrs Aj and Wootters's basis can transforms under 5*0(4, c) into a genèric 2x2 
matrix with the same set of positive numbers but with new Wootters's basis, where the local 
unitary transformations correspond to 50(4, r) transformations, hence, p can be represented 
as coset space 50(4, c)/50(4, r) together with positive numbers Aj. By giving an cxplicit 
parameterization we characterize a genèric orbit of group of local unitary transformations. 

Keywords: Quantum entanglement, Concurrence , Lewenstein-Sanpera decom- 
position, Two qubit systems, Wootters's basis, 50(4, c) 
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1 Introduction 

Perhaps, quantum entanglement is the most non-classical features of quantum mechanics [1, 2] 
which has recently been attracted much attention although it was discovered many dècades ago by 
Einstein and Schròdinger [1, 2]. It plays a central role in quantum information theory and provides 
potential resource for quantum communication and information processing [3, 4, 5]. Entanglement 
is usually arise from quantum correlations between separated subsystems which can not be created 
by local actions on each subsystems. By definition, a bipartite mixed state p is said to be separable 
if it can be expressed as 



where p\ and p\ denote density matrices of subsystems 1 and 2, respectively. Otherwise the 
state is entangled. 

The central tasks of quantum information theory is to characterize and quantify entangled 
states. A first attempt in characterization of entangled states has been made by Peres and Horodecki 
family [6, 7]. Peres showed that a necessary condition for separability of a two partite system is 
that its partial transposition be positive. Horodeckis have shown that this condition is sufficient 
for separability of composite systems only for dimensions 2 (g> 2 and 2 <g) 3. 

There is also an increasing attention in quantifying entanglement, particular ly for mixed states 
of a bipartite system, and a number of measures have been proposed [5, 8, 9, 10]. Among them the 
entanglement of formation has more importance, since it intends to quantify the resources needed 
to create a given entangled state. 

An interesting description of entanglement is Lewenstein-Sanpera decomposition [11]. Lewen- 
stein and Sanpera in [11] showed that any two partite density matrix can be represented optimally 



P = w i Pi ® Pi 



Wi>0, ^2 Wi = 1 
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as a sum of a separable state and an entangled state. They have also shown that for two qubit 
systems the decomposition reduces to a mixture of a mixed separable state and an entangled pure 
state, thus all non-separability content of the state is concentrated in the pure entangled state. 
This leads to an unambiguous measure of entanglement for any two qubit state as entanglement of 
pure state multiplied by the weight of pure part in the decomposition. 

The numerical method for finding the BSA has been reported in Ref. [11]. Also in two qubit 
systems some analytical results for special states were found in [12]. An attempt to generalize the 
results of Ref. [11] is made in [13]. In [14] an algebraic approach to find BSA of a two qubit state 
is attempted. They have also shown that the weight of the entangled part in the decomposition 
is equal to the concurrence of the state. In [15] we have obtained an analytical expression for L-S 
decomposition of Bell decomposablc (BD) states. We have also obtained the optimal decomposition 
for a class of states obtained from BD states via some LQCC actions. 

Wootters in [10] has shown that for any two qubit density matrix p there always exist a de- 
composition p = J2í\ x í) ( x i\ such that (xi\xj) = XiSij, where Àj are square roots of eigenvalues, 
in decreasing order, of the non-Hermitian matrix pp. Based on this the concurrence of the mixed 
state p is defined by max(0, Ai — À2 — A3 — À4) [10]. 

In this paper, by using Wootters's basis, we give an analytical expression for optimal Lewenstein- 
Sanpera decomposition for any two qubit density matrix. We show that all entanglement content 
of the state is concentrated in the Wootters's state \x\) associated with the largest eigenvalue Ai. 
It is also shown that the average concurrence of the decomposition is equal to the concurrence of 
the state. 

It is shown that a given density matrix p with corresponding set of positive numbers Aj and 
Wootters's basis can transforms under SO(4, c) into a genèric 2x2 matrix with the same set of 
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positive numbers but with new Wootters's basis, where the local unitary transformations correspond 
to ,SO(4,r) transformations, hence, p can be represented as coset space £0(4, c) /SO (4, r) together 
with positive numbers Aj. By giving an explicit parameterization we characterize a genèric orbit of 
group of local unitary transformations. 

The paper is organized as follows. In section 2 we review concurrence for two qubit density 
matrix. In section 3 we give an analytical expression for L-S decomposition of a genèric density 
matrix. We prové that the decomposition is optimal. Characterization of the density matrix in 
terms of orthogonal group is presented in section 4. Section 4 is devoted to explicit parameterization 
of a genèric density matrix up to a local unitary transformation. The paper is ended with a brief 
conclusion in section 5. 

2 Concurrence 

In this section we first review concurrence of two qubit mixed states. From the various measures 
proposed to quantify entanglement, the entanglement of formation has a special position which 
in fact intends to quantify the resources needed to create a given entangled state [5]. In the case 
of pure state if the density matrix obtained from partial trace over other subsystems is not pure 
the state is entangled. For the pure state of a bipartite system, entropy of the density matrix 
associated with either of the the two subsystems is a good measure of entanglement 

E(ip) = —Tr(pA log2 Pa) = -Tr(p B log 2 p B ), 

where p\ = 2>b(|V>) (V'l) an d Pb defined similar ly. Due to classical correlations where exist in the 
mixed state each subsystem can have non-zero entropy even if there is no entanglement, therefore 
von Neumann entropy of a subsystem is no longer a good measure of entanglement. For a mixed 
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state entanglement of formation is defined as the minimum of average entropy of the state over all 
pure state decompositions of the state 

E f (p) = mmJ2piE^i). (2-1) 

i 

Wootters in [10] has shown that for a two qubit system entanglement of formation of a mixed 
state p can be defined as 

E f {p) = H (± + \y/ï=&) , (2-2) 
where H(x) = — xlnx — (1 — x) ln (1 — x) is binary entropy and concurrence C{p) is defined by 

C(p) = max{0, Ai - A 2 - A 3 - A 4 }, (2-3) 

where the \ are the non-negative eigenvalues, in decreasing order, of the Hermitian matrix R = 



yfPPyfP and 

p = (a y <g> (Ty)p*((Ty ® (Ty), (2-4) 

where p* is the complex conj ugate of p when it is expressed in a Standard basis such as{|TT))|íl)}>{|lí))|ll)} 
and a y represent Pauli matrix in local basis {|T) > II)} • 

Consider a genèric two qubit density matrix p with its subnormalized orthogonal eigenvectors 
\ví), i.e. p = J2i \ v i) ( v i\- There always exist a decomposition [10] 

P = ^\xí)(xí\ (2-5) 

i 

where Wootters states |xj) are defined by 

4 

I**) = E U tj h) . for » = 1,2,3,4, (2-6) 

j 

such that 

( Xi | Xj) = {UrU T ) tj = XiSij, (2-7) 
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where = (ví \ Vj) is a symmetric but not necessarily Hermitian matrix. To construct \xí) we use 
the fact that for any symmetric matrix r one can always find a unitary matrix U in such a way that 
Xi are real and non-negative, that is, they are the square roots of eigenvalues of rr* which are the 
same as eigenvalues of R. Alternatively one can always find U such that Aj appear in decreasing 
order. 

3 Lewenstein-Sanpera decomposition 

According to Lewenstein-Sanpera decomposition [11], any two qubit density matrix p can be written 

as 



where p sep is a separable density matrix and is a pure entangled state. The Lewenstein-Sanpera 
decomposition of a given density matrix p is not unique and, in general, there is a continuum set 
of L-S decomposition to choose from. The optimal decomposition is, however, unique for which A 
is maximal, and 



In Ref. [11] Lewenstein and Sanpera have shown that any other decomposition of the form p = 
Xp se p + (1 - À)|^)(^| , with À G [0, 1] such that p / pM) necessarily implies that À < \^ [11]. 

Here in this section we obtain L-S decomposition for a genèric two qubit density matrix by 
using Wootters states. First we define states |a^) in terms of Wootters states as 




A e [o, i] 



(3-8) 





(3-9) 




Xi) 



for i = 1,2,3,4. 



(3-10) 



v/A? 
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Thus the decomposition given in Eq. (2-5) becomes 

p^a.KX^i, (3-n) 

i 

where can be rewritten in the following form 



(Ai - A 2 - À 3 - A 4 ) {x[\ + (A 2 + A 3 + A 4 ) \x[) {x[\ + El 2 ^ 



X > )(x> 



(3-12) 



;j=2 "3 
= (1-A)|V) {lp\+\p S e P , 

where separable density matrix p sep and entangled pure state \ip) are given by 

Psep = + \\* ^ ) \ Xl "> ( Xl \ + \ S l x J> fol ' ( 3 " 13 ) 

and 

IY>> = "7==, (3-14) 



respectively, and parameter A is equal to 



A = l-(' Al - A2 " A3 ~ V ) ( Xl \ Xl ). (3-15) 



Ai 

Equation (3-14) shows that all entanglement content of the state is concentrated in the Wootters 
state \xi) associated to the largest eigenvalues Ai. It is also worth to note that average concurrence 
of the decomposition is 

(1 - A) = (Ai - A 2 - A 3 - A 4 ), (3-16) 

that is, it is equal to the concurrence of the state. 

In order to show that the decomposition given in (3-12) is optimal we first show that p sep can 
be represented as a convex sum of product states. First note that p sep can be written as 

4 

/Vp = £K><*íI> ( 3 - 17 ) 
i=l 



L-S decomposition for 2x2 systems 
where \x'() are defined by 



A2 + A3 + À4 



1 



ï) = \/^-j^ 4> = — forj= 2,3,4. (3-18) 



Va 



Obviously, the basis x") (i = 1,2,3,4) satisfy the following relations 

where (j = 1, 2, 3, 4) are corresponding dual basis. 

Wootters in [10] has shown that any two qubit separable density matrix with decomposition 
given in Eq. (3-17) can be expanded in terms of following product states 



\x'{)+é^ 


\4) + e iea 


|4) + e m 


14)) , 


(3-19) 


\ x '[) + é 6 * 


|4) - e Wa 


|4) - e ie * 


14)) , 


(3-20) 


\x'[) - e id2 


\x'ï) + é 9 * 


|4) - 


14)) > 


(3-21) 


\x'[) - e ie > 




|4) + é e ± 


14)) . 


(3-22) 



1 z, 2 v 

w - è (■ 

N> = j (' 

where zero concurrence is guaranteed with J2j=i e 2l ^X" = 0. 

Now using the fact that the eigenvalues \" of boundary separable states p sep satisfy constraint 
X'l — X'2 — A3 — À4 = 0, we can choose phase factors Oi as 62 = #3 = 64 = 0\ + tt. Choosing 0\ = 
we arrive at the following product ensemble for p sep 



^ = 27x (y 2 + Ai + ^ ^ ~ i ' X2 ^ + * ' X3 ^ + i ' X4 w ' ( - 3 " 24 ' ) 



= 27X (y 2 + Ai + A4 ^ ~ i \ x 3) + i \ x 'òj , (3-25) 
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= ( V A2 + A 3 + ^ 4 + * + * ' X3 ^ ~~ Í ' X4 ^ j ' (3-26) 
It can be easily seen that all \zí) have zero concurrence and also p sep can be expanded as 

4 

Psep = ^2\zí) (Zi\ . (3-27) 

1=1 

In the rest of this section we will prové that the decomposition (3-12) is the optimal one. To 
do so we have to find a decomposition for p sep in terms of product states \e a , f a ), i.e. 

Psep = ^ ' A Q \&cti fa) (fia; Ía| (3-28) 
a 

such that the following conditions are satisfied [11] 

i) All A a are maximal with respect to p a = A a \e a , f a ) (e a , f a \ + (1 — A) and projector 

Pa = | &a i fa) (Cq , /a | • 

ii) All pairs (A Q , A^) are maximal with respect to p a p = A a \e a , f a ) {e a , f a \+Ap \ep, fp) (ep, fp\ + 
(1 — À) \tp) (ip\ and the pairs of projector (P a , Pp). 

Then according to [11] p sep is BSA and the decomposition given in Eq. (3-12) is optimal. 
Lewenstein and Sanpera in [11] have shown that A a is maximal with respect to p a and P a = 
\e a ,f a ) (tp\ iff 

a) if \e a , f a ) G" TZ(pa) then A a = 0, and b) if \e a , f a ) G TZ(p a ) then A a = (e a , f a \ p~ l \e a , fa}' 1 > 
0. They have also shown that a pair (Ai, A2) is maximal with respect to pu and a pair of projectors 
(Pi,P 2 ) iff: 

a) if \e\,fi), |e2,/2) do not belong to IZ(pn) then Ai = A2 = 0; b) if \ei,f\) does not belong, 
while |e 2 ,/ 2 ) G T^(pn) then A x = 0, A 2 = (e 2 , f 2 \ PÏ2 1^2, /2>~ 1 ; c) if |ei,/i), |e 2 ,/ 2 ) G TZ(pn) and 
(ei, /i| l e 2, ^2) = then Aj = (e*, fc\ p^ |e Í5 /i) -1 , i = 1, 2; d) finally, if |ei, /1) , |e 2 , / 2 ) G 7£(pi 2 ) 
and (ei, /i| p^ 1 \^2, h) ^ ^ then 



L-S decomposition for 2x2 systems 



11 



Ai = ((e 2 , f 2 \ P12 |e 2 , Í2) - | (ei, /i| p 12 X |e 2 , / 2 ) |)/£>, 

(3-29) 

A 2 = «ei, /i| p^ 1 l e i> /1) - I (ei, /i| P^ 1 l e 2, /2> \)/D, 
where L> = (ei,/i| p^ 1 |ei,/i) (e 2 ,/ 2 | p^ 1 |e 2 ,/ 2 ) - | (ei,/i| p_i |e 2 ,/ 2 ) | 2 . 

Now we return to show that the decomposition given in (3-12) is optimal. We first consider the 
cases that p has full rank. Let us consider the set of four product vectors {\zí)} and one entangled 
state \x\). In Ref. [10] it is shown that the ensemble {\zí}} are linearly independent, also it is 
straightforward to see that three vectors \z a ) , \zp) and \x\) are linearly independent. Now let us 
consider matrices p a = A a \z a ) (z a \ + (1 — À) \x±) {x±\. Due to independency of \z a ) and \x\) we can 
deduce that the range of p a is two dimensional, thus after restriction to its range and defining dual 
basis \z a ) and \x\) we can expand restricted inverse p" 1 as p" 1 = A^ 1 )^ } ( z a \ + (1 — À) _1 |fi }(x\\ 
(see appendix). Using Eq. (5-59) it is easy to see that {z a \p^ t l \z a ) = A~ l . This shows that A Q are 
maximal with respect to p a and the projector P a = \z a ) (z a \. 

Similarly considering matrices p a p = A a \z a ) (z a \ + \zp) (z@\ + (1 — À) \x±) (xi\ and considering 
the independency of vectors \z a ), \zp) and \x±) we see that the rang of p a p is three dimensional where 
after restriction to its range and defining their dual basis \z a ), \zp) and we can write restricted 
inverse p~p as p~p = A~ 1 \z a )(z a \ + A^\zp )( zp\ + (1 — A) _1 |xi )( x\\. Then it is straightforward 
to get (z a \ p~p \z a ) = A" 1 , (zp\ p"^ \zp) = A^ 1 and (z a \ p~^ \zp) = 0. This implies that the pair 
(A a ,Ap) are maximal with respect to p a p and the pair of projectors (P Q ,P^), thus complete the 
proof that the decomposition given in Eq. (3-12) is optimal. 

We now consider the cases that p has rank three, that is A4 = 0. In this case the pairs 
{\ z i) 1 \ z i)} and also {|z 2 ) > | ^3) } are no longer independent with respect to \x\). In former case 
we can evaluate \x\) in terms of \z\) and \z±) then matrix P14 can be written in terms of two 
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basis \zi) and \za) which yields after some calculations, ^zi|/9 14 1 |zi^ = ^ (A4 + (1 — À) ( a 2 +a 3 
(z4\pÏÏ\zí) = £ (Ai + (1 - A) (^)) and (z^W) = j£ ((1 - A) (^)) where T 14 = 
(AiA 4 + (Ai + A 2 )(l - A) ( a 2 +a 3 ))- Using the above results together with Eqs. (3-29) we ob- 
tain the maximality of pair (A^A^ with respect to p±4 and the pair of projectors (Pi,P 4 ). 
Similar ly in the second case one can express \x\) in terms of \z 2 ) and \z 3 ) and evaluate p 23 , 
which get (z 2 \ P2 i\z 2 ) = r^(A3 + (l-A)(^)), {z 3 \p 23 l \z 3 ) = £ (a 2 + (1 - A) (^)) 
and (z 2 \p^\z 3 ) = r i((l_A)(^)) with T 23 = (a 2 A 3 + (A 2 + A 3 )(l - A) , to- 

gether with Eqs. (3-29) we obtain the maximality of pair (A 2 ,A 3 ) with respect to p 23 and the pair 
of projectors (P 2 , P 3 ). For other choices of a and (3 three vectors \z a ) ,\zp) and \x\) remain linearly 
independent thus we can prové maximality of pairs (A Q , A^) in the same way that we proved it in 
full rank case. 

Finally let us consider cases that p has rank two, that is A 3 = A4 = 0. In this case we have 
\zi) = \z 2 ) and \z 3 ) = \z 4 ). It is now sufficient to take \z\) and \z 3 ) as product ensemble. But in 
this case vectors \z\) and \z 3 ) are not independent with respect to \xi). We express \x±) in terms 
of \zi) and \z 3 ) then matrix p\ 3 can be written in terms of two vectors \z\) and \z 3 ) and after some 
calculations we get ( Zl \p^\ Zl ) = ^ (a 3 + (1 - A) (^)), (z 3 \p^\z 3 ) = ^ (a 4 + (1 - A) (^)) 
and ( Zl \p^\z 3 ) = ^ ((1 - A) (£)) where T 13 = (A1A3 + (Ai + A 3 )(l - A) (£)). Using the 
above results together with Eqs. (3-29) we obtain the maximality of pairs (Ai, A 3 ) with respect to 
pi 3 and the pairs of projectors (P\,P 3 ). 
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4 Coset structure for a genèric 2x2 density matrix in Wootters's 
basis 

In this section we obtain an explícit parameterization for a genèric two qubit density matrix in 
Wootters's basis. To this aim for any density matrix p with decomposition given in Eq. (3-11) we 
define matrix X as 

X=(\x' 1 ),\x' 2 ),\x' 3 ),\x' 4 )). 



Analogously by defining matrix 

X = [ x'i) , x'i) , x' 3 ) , x' 4j j . 
Eq. (2-7) takes the following form 



(4-30) 



(4-31) 



X*X = X T a y ® G y X = I. 



(4-32) 



Since matrix a y <S> <r y is symmetric it can be diagonalized as 



a y ®(jy = T rj 2 0, 



(4-33) 



where O is an orthogonal matrix defined by 



V2 



1 








1 





1 


1 








1 


-1 





1 








-1 



(4-34) 
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1 











i 













as 







(4-35) 



Y T Y = L 



(4-36) 



where 



Y = r] OX. (4-37) 

Equation (4-36) shows that Y is a complex 4-dimensional orthogonal matrix. This means that a 
given density matrix p with corresponding set of positive numbers Aj and Wootters's basis can 
transforms under SO '(4, c) into a genèric 2x2 density matrix with the same set of positive numbers 
but with new Wootters's basis. This implies that the space of two qubit density matrices can be 
characterize with 12-dimensional (as real manifold) space of complex orthogonal group 50(4, c) 
together with four positive numbers Aj. Of course the normalization condition reduces number of 
parameters to 15. 

As far as entanglement is concerned the states p and p' are equivalent if they are on the same 
orbit of the group of local transformation, that is, if there exist local unitary transformation U\®Ui 
such that p' = {Ui®U2)p(Ui®U2)\ where U\ and Ui are unitary transformations acting on Hilbert 
spaces of partides A and B, respectively. 

It can be easily seen that under above mentioned local unitary transformations of density matrix 
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p, the matrix X transforms as 

X -> X' = (t/i ® C/ 2 )X (4-38) 

It is worth to mention that X' also satisfy Eq. (4-32). To show that this is indeed the case, we 
need to note that X' T u y <8> a y X' = X T (Uf a y Ui) <8> {U^ cr y U2)X . By using (a y )ij = —iíij we get 
{U T a y U)ij = —iekiUkiUij = — idet(U)eij = a y , where the fact that Ui G 517(2), thus having unit 
determinant, have been used. This implies that 

X' T a y ®a y X' = I. (4-39) 

By defining Y' as 

Y' = t\OX', (4-40) 
one can easily show that Y' is also satisfied orthogonality condition 

Y' T Y' = I. (4-41) 

Now by using Eq. (4-40) and inverting Eq. (4-38), we can express Y' in terms of Y 

Y' = (t/0)(Z7i ^2)^0)-^. (4-42) 

Now by using the fact that (rjO) exp(Ui (3lA2){riO)~ l = exp((rj 0){U\ ® U-i){j] 0)~ v ) and using 
the explicit form for generators (U\ <g> W2) of local group, one can after some algebraic calculations 
see that (rjO)(Ui (8> 1^2) {v 0)~ l is real antisymmetric matrix. This means that under local unitary 
transformations matrix Y transforms with 50(4, r) group. So we can parameterize the space of 
two qubit density matrices as 6-dimensional coset space 50(4, c)/50(4, r) together with 4 positive 
numbers Àj, which again normailzation condition reduces the number of parameters to 9. 
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In the following we will obtain an explícit parameterization for a genèric two qubit density 
matrix. First note that we can decompose coset SO(4,c)/SO(4,r) as [16] 



SO (4, c) 



SO(4,c)/SO(4,r) 



fSO(2,c) SO(2,c) 



(4-43) 



SO(4,r) SO(2,c)/SO(2,r) SO{2,c)/SO{2,r) ^ \SO{2,r) SO(2,r)J ' 

that is, coset representative Y can be decomposed as Y = Y{Y2- One can easily show that coset 
representative of SO(2,c)/SO(2,r) has the following form 

/ 



exp 



V 



o i4> 

-icfi 



\ ( \ 

coshò ismhò 



) 



V 



-i sinh (fi cosh (fi 



(4-44) 



Thus Y2 can be written as 



Yo 



On the other hand Y\ can be evaluated as 



cosh <fii i sinh <fii 









—i sinh (fií cosh (fií 






cosh (fi>2 i sinh (fi>2 









—i sinh 4>2 cosh 02 



(4-45) 



Yi = exp 






iB 


-ÍB T 






cosh VBB T 



„• sinh V B T B j^T 



sinh V B T B 
\fWB 



cosh Vb^B 





1 







(4-46) 



where i? is a 2 x 2 matrix and in the last step we used C = B smh , ^J 5 ^ . Now using the singular 

y B B 



value decomposition C = O1DO2 ', Eq. (4-46) becomes 



OiVl + D 2 Oj 



V 



-io 2 Dó{ 



iOiDOl 



O2VTTDÏOÏ 



(4-47) 



where D is a non-negative diagonal matrix. It can be easily seen that Eq. (4-47) can be decomposed 
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as 



Y 1 



Oi 








o 2 



( 



I 



V 



-íd 



íd 



VI + D 2 



By combining Eqs. (4-45) and (4-48) we get 



Y 



( 


°1 


/ 






Oi 


V/ + D 2 


íd 


{ 


o 2 J 


V 


-ÍD 


Vl + D 2 



oï 








ol 










o 2 



Finally using parameterization given in Eq. (4-45) we get 



Y = 



(4-48) 



(4-49) 



/ 


cosh Oi i sinh Q\ 





\ 


/ 


cosh Çi 





i sinh 







—i sinh 6\ cosh Q\ 









cosh Ç2 





i sinh Ç2 







cosh 02 i sinh 6*2 






— isinh £i 





cosh £1 





V 


— i sinh # 2 cosh (9 2 


/ 


V 





-i sin Ç 2 





cosh Ç2 



cosh <pi i sinh 0i 









—i sinh 0i cosh 0i 






cosh 02 i sinh 02 









—i sinh 02 cosh 02 



(4-50) 



where sinhÇj (for í = 1,2) are diagonal elements of D with the conditions & > 0. 
Using above results and Eq. (4-30) and (4-37) we can evaluate the states \xí) as 



Fi 



— (sinh£i sinh #2 cosh0i + sinh £2 coshé>2 sinh0i) — í(cosh£i coshé*i cosh0i + cosh £2 sinhé>i sinh0i) 

— (sinh£i cosh # 2 cosh0i + sinh£ 2 sinh# 2 sinh0i) — i(cosh£i sinh#i cosh0i + cosh Ç 2 cosh 0\ sinh0i) 

(sinh £ 1 cosh 62 cosh 0i + sinh £2 sinh 62 sinh 0i ) — i (cosh £1 sinh 9\ cosh 0i + cosh £2 cosh 9\ sinh </>i ) 

(sinh £1 sinh 62 cosh 0i + sinh £ 2 cosh 62 sinh 0i ) — i (cosh £1 cosh 81 cosh 0i + cosh £ 2 sinh 0\ sinh 0i ) 

(4-51) 
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1^2) = 



I 2 *} = \ hr 



( \ 

(cosh £1 cosh Q\ sinh 0i + cosh £2 sinh 6\ cosh (f>i) — i (sinh £1 sinh 02 sinh 0i + sinh £2 cosh 02 cosh 0i ) 

(cosh £1 sinh 0i sinh 0i + cosh £2 cosh d\ cosh 0i ) — i (sinh £1 cosh #2 sinh 0i + sinh £2 sinh 02 cosh 0i ) 

(cosh£i sinh 0i sinh0i + cosh£ 2 cosh0i cosh0i) + i(sinh£i cosh0 2 sinh0i + sinh£ 2 sinh0 2 cosh0i) 

(cosh£i cosh0i sinh0i + cosh£ 2 sinh0i cosh0i) + i(sinh£i sinh0 2 sinh0i + sinh£ 2 cosh0 2 cosh0i) 

(4-52) 

(sinh £1 cosh 0i cosh 2 + sinh £ 2 sinh 0i sinh 2 ) — i (cosh £1 sinh 62 cosh 2 + cosh £ 2 cosh 2 sinh 02 ) 
(sinh£i sinh 0i cosh 02 + sinh £2 cosh 9\ sinh 02) — ï(cosh£i cosh 02 cosh 02 + cosh £2 sinh 02 sinh 02) 
(sinh £1 sinh 0i cosh 2 + sinh £ 2 cosh 0i sinh 2 ) + i (cosh £1 cosh 2 cosh 2 + cosh £ 2 sinh 2 sinh 2 ) 



/ 



(sinh £1 cosh 0i cosh 2 + sinh £ 2 sinh 0i sinh 2 ) + i (cosh £1 sinh 2 cosh 2 + cosh £ 2 cosh 2 sinh 2 ) 

(4-53) 

(coshÇi sinh 2 sinh0 2 + cosh£ 2 cosh0 2 cosh0 2 ) + ï(sinh£i cosh0i sinh0 2 + sinh£ 2 sinh0i cosh 02) 
(cosh^i cosh 2 sinh 02 + cosh£ 2 sinh0 2 cosh0 2 ) + i(sinh£i sinh0i sinh0 2 + sinh£ 2 cosh0i cosh0 2 ) 
— (cosh Çi cosh 02 sinh 02 + cosh £2 sinh 02 cosh 02 ) + i (sinh £1 sinh 0i sinh 02 + sinh £2 cosh 0i cosh 02 ) 



— (cosh £1 sinh 2 sinh 2 + cosh £ 2 cosh 2 cosh 2 ) + i (sinh £ 1 cosh 0i sinh 2 + sinh Ç 2 sinh 0i cosh 2 ) 

(4-54) 



Equations (4-51) to (4-54) together with normalization condition Ylí=i i x i\ x i) = 1 gi ye a parame- 
terization for a genèric orbit of two qubit density matrix up to local unitary group. As an examplc 
let us consider Bell decomposable states p = J2í=iPí IV'j) (V'íIj where Pi > 0, J2iPi = 1- F° r these 
states by choosing 9\ = 62 = £1 = Ç2 = 01 = 4>2 = we get Àj = pi and states \xí) are given by 



Xl) = -ïy 


/pïlV'i) > 


IV>1> = 


^(ITT) + IU)), 


(4-55) 


M = A 


/P2 |^2) , 


|V<2> = 


^(ITD + IIT)), 


(4-56) 


^3) = ~i\ 


/P3 1^3) , 


IV>3> = 


^(ITD-IIT)), 


(4-57) 


|x 4 ) = 


s/pi |V>4> , 


|Y>4> = 


^(ITT)-IU)). 


(4-58) 
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5 Conclusion 

We have obtained Lewenstein-Sanpera decomposition for a genèric two qubit density matrix by 
using Wootters's basis. It is shown that the average concurrence of the decomposition is equal 
to the concurrence of the state. It is also shown that all entanglement content of the state is 
concentrated in the Wootters's state \x{) associated with the largest eigenvalue Ai. It is shown 
that a given density matrix p with corresponding set of positive numbers Aj and Wootters's basis 
can transforms under SO(4, c) into a genèric 2x2 matrix with the same set of positive numbers but 
with new Wootters's basis. We have also shown that the local unitary transformations correspond 
to £0(4, r) transformations, hence, p can be represented as coset space SO(4,c)/SO(4,r) together 
with positive numbers Aj. 
Appendix 

Let us consider the set of linear ly independent vectors {]&)}, then one can define their dual 





(5-59) 



hold. It is straightforward 




following 



completeness relation 



Y,\k){4>i\ = i, J2\4> i )(è i \ = i. 



(5-60) 



Consider an invertible operator M which is expanded in terms of states \4>i) as 




(5-61) 



i 



Then the inverse of M denoted by M 



can be expanded in terms of dual bases as 



M 



i 



^2b ij \4> i )(4> j 



(5-62) 
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where bij = (A^ 1 )^ and A; L j = aij. 
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